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Abstract. In this paper we present the applications of methods from wavelet analysis 
to polynomial approximations for a number of accelerator physics problems. In the 
general case we have the solution as a multiresolution expansion in the base of com- 
pactly supported wavelet basis. The solution is parametrized by the solutions of two 



, 

reduced algebraical problems, one is nonlinear and the second is some linear problem, 
i <~| which is obtained from one of the next wavelet constructions: Fast Wavelet Transform, 

Stationary Subdivision Schemes, the method of Connection Coefficients. According to 
the orbit method and by using construction from the geometric quantization theory we 
construct the symplectic and Poisson structures associated with generalized wavelets 
!> ■ by using metaplectic structure. We consider wavelet approach to the calculations of 

Melnikov functions in the theory of homoclinic chaos in perturbed Hamiltonian systems 
and for parametrization of Arnold- Weinstein curves in Floer variational approach. 

o 

S> ; INTRODUCTION. 

o 



In this paper we consider the following problems: the calculation of orbital mo- 
tion in storage rings, some aspects of symplectic invariant approach to wavelet 
computations, Melnikov functions approach in the theory of homoclinic chaos, the 
calculation of Arnold- Weinstein curves (periodic loops) in Hamiltonian systems. 
The key point in the solution of these problems is the use of the methods of wavelet 
analysis, relatively novel set of mathematical methods, which gives us a possibility 
to work with well-localized bases in functional spaces and with the general type 
of operators (including pseudodifferential) in such bases. Our problem as many 
related problems in the framework of our type of approximations of complicated 
physical nonlinearities is reduced to the problem of the solving of the systems of dif- 
ferential equations with polynomial nonlinearities with or without some constraints. 
In this paper we consider as the main example the particle motion in storage rings 
in standard approach, which is based on consideration in [1], [2]. Starting from 
Hamiltonian, which described classical dynamics in storage rings 



H{r,P,t) = c {n 2 + m 2 c 2 } 1/2 + e 



and using Serret-Frenet parametrization, we have the following Hamiltonian for 
orbital motion in machine coordinates: 

H(x,p x ,z,p z ,a,p a ; s) = p a - [1 + f{p a ] ■ [1 + K x ■ x + K z ■ z] x (2) 
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Then, after standard manipulations with truncation of power series expansion 
of square root we arrive to the following approximated Hamiltonian for particle 
motion: 
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and the corresponding equations of motion: 
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Then we use series expansion of function f{p a ) from [2]: 



f( P ,) = f(0) + f'(0)p a + f"(0)-pl + 
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and the corresponding expansion of RHS of equations (4). In the following we 
take into account only an arbitrary polynomial (in terms of dynamical variables) 
expressions and neglecting all nonpolynomial types of expressions, i.e. we consider 
such approximations of RHS, which are not more than polynomial functions in 
dynamical variables and arbitrary functions of independent variable s ("time" in 
our case, if we consider our system of equations as dynamical problem). 



I POLYNOMIAL DYNAMICS 
Introduction. 

The first main part of our consideration is some variational approach to this 
problem, which reduces initial problem to the problem of solution of functional 
equations at the first stage and some algebraical problems at the second stage. We 
consider also two private cases of our general construction. In the first case (par- 
ticular) we have for Riccati equations (particular quadratic approximations) the 
solution as a series on shifted Legendre polynomials, which is parameterized by the 
solution of reduced algebraical (also Riccati) system of equations. This is only an 
example of general construction. In the second case (general polynomial system) 
we have the solution in a compactly supported wavelet basis. Multiresolution ex- 
pansion is the second main part of our construction. The solution is parameterized 
by solutions of two reduced algebraical problems, one as in the first case and the 
second is some linear problem, which is obtained from one of the next wavelet con- 
struction: Fast Wavelet Transform (FWT), Stationary Subdivision Schemes (SSS), 
the method of Connection Coefficients (CC). 



Variational method. 



Our problems may be formulated as the systems of ordinary differential equations 
dxi/dt = fi(xj,t), = 1, ...,n) with fixed initial conditions Xi(0), where fi are 



not more than polynomial functions of dynamical variables Xj and have arbitrary 
dependence of time. Because of time dilation we can consider only next time 
interval: < t < 1. Let us consider a set of functions = Xidi/i/dt + fiUi 

and a set of functional F^x) = J$ &i(t)dt — x^i |q, where yi(t)(yi(0) = 0) are 
dual variables. It is obvious that the initial system and the system Fi(x) = are 
equivalent. In the last part we consider the symplectization of this approach. Now 
we consider formal expansions for Xj,^: 

X l (t)=X l (0)+Y, X *Mt) y J {t)=T.V>r(t) ) (5) 

k r 

where because of initial conditions we need only (pk(0) = 0. Then we have the 
following reduced algebraical system of equations on the set of unknown coefficients 
Af of expansions (5): 

^^A fc -7[(A,) = (6) 
k 

Its coefficients are fi^r = Jo Pkitf^ritydt, 7[ = Jq 1 fi(xj,t)(p r (t)dt. Now, when we 
solve system (6) and determine unknown coefficients from formal expansion (5) 
we therefore obtain the solution of our initial problem. It should be noted if we 
consider only truncated expansion (5) with N terms then we have from (6) the 
system of N x n algebraical equations and the degree of this algebraical system 
coincides with degree of initial differential system. So, we have the solution of the 
initial nonlinear (polynomial) problem in the form 

N 

x t (t)=x t (0) + Y, X *Mt), (7) 

k=l 

where coefficients are roots of the corresponding reduced algebraical problem (6). 
Consequently, we have a parametrization of solution of initial problem by solution 
of reduced algebraical problem (6). But in general case, when the problem of com- 
putations of coefficients of reduced algebraical system (6) is not solved explicitly as 
in the quadratic case, which we shall consider below, we have also parametrization 
of solution (4) by solution of corresponding problems, which appear when we need 
to calculate coefficients of (6). As we shall see, these problems may be explicitly 
solved in wavelet approach. 

The solutions 

Next we consider the construction of explicit time solution for our problem. 
The obtained solutions are given in the form (7), where in our first case we have 
Xk(t) = Qk(t), where Qk{t) are shifted Legendre polynomials and A^. are roots 
of reduced quadratic system of equations. In wavelet case Xk(t) correspond to 



multiresolution expansions in the base of compactly supported wavelets and X l k are 
the roots of corresponding general polynomial system (6) with coefficients, which 
are given by FWT, SSS or CC constructions. According to the variational method 
to give the reduction from differential to algebraical system of equations we need 
compute the objects 7^ and /iji, which are constructed from objects: 



J X>(T)X j (T)dT = a j F 1 (i : 0) + F 1 (i,j), (8) 
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if j + k + l = 2m,me Z,and a klj = if j + k + l = 2m+l; where R(i) = (2i)\/(2H\) 2 , 
Qi = &i + P*, where the second equality in the formulae for a,u, /i, f3,a hold for 
the first case. 



Wavelet computations. 

Now we give construction for computations of objects (8) in the wavelet case. 
We use some constructions from multiresolution analysis: a sequence of successive 
approximation closed subspaces Vf ...V 2 C V\ C Vo C V_i C V_ 2 C ... satisfying 

the following properties: f| Vj = 0, [jVj = L 2 (R), /(x) G Vj <=> f(2x) E V j+1 

jez jez 

There is a function (p e Vo such that {vo,*;^) = ~~ ^)fcgz) f° rms a Riesz basis 
for Vq. We use compactly supported wavelet basis: orthonormal basis for functions 
in L 2 (R). As usually (p(x) is a scaling function, ip(x) is a wavelet function, where 
(Pi(x) = ip(x — i). Scaling relation that defines ip, if) are 

N-l N-l N-2 

f( x ) = H a k <p(2x - k) = a k ip k (2x), ip(x) = ^2 *'a k+1 ip(2x + k) 

k=0 k=0 k=-l 

Let be / : R — > C and the wavelet expansion is 
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f( x ) = S c We( x ) + H c jkipjk(x) (9) 
tez j=ofcez 



The indices k, £ and j represent translation and scaling, respectively 



tpflix) = 2^V(2^x - e),4 jk (x) = 2i/ 2 iP(Vx - k) 

The set {(pj t k}kez forms a Riesz basis for Vj. Let Wj be the orthonormal complement 
of Vj with respect to Vj + \. Just as Vj is spanned by dilation and translations of 
the scaling function, so are Wj spanned by translations and dilation of the mother 
wavelet ipjk(x). If in formulae (9) Cjk = for j > J, then f(x) has an alternative 
expansion in terms of dilated scaling functions only f(x) — cje(pje(x). This is 

a finite wavelet expansion, it can be written solely in terms of translated scaling 
functions. We use wavelet ip(x), which has k vanishing moments / x k ip(x)d(x) = 0, 
or equivalently x k = J2ce(pe(x) for each k, < k < K. Also we have the shortest 
possible support: scaling function DN (where A" is even integer) will have support 
[0, N — 1] and N/2 vanishing moments. There exists A > such that DN has AA" 
continuous derivatives; for small N,\> 0.55. To solve our second associated linear 
problem we need to evaluate derivatives of f(x) in terms of (fi(x). Let be ip™ = 
d n Lp i{x) j dx n . We derive the wavelet - Galerkin approximation of a differentiated 
f(x) as f d (x) = J^eCWiix) an d values (pg(x) can be expanded in terms of (p(x) 

oo 

= ^2^m(Pm(x), A m = / (p d (x)(p m (x)dx 
m 

The coefficients X m are 2-term connection coefficients. In general we need to find 
(di > 0) 

oo 

J II ( 10 ) 
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For Riccati case we need to evaluate two and three connection coefficients 

oo 

Af d2 = J% d '(x)^(x)dx, A d ^ = j <ffl(x)<pt>(x)<f%(x)dx 

— oo 

According to CC method [12] we use the next construction. When A" in scaling 
equation is a finite even positive integer the function (p(x) has compact support 
contained in [0, N — 1]. For a fixed triple (di, d 2 , ^3) only some A d ^ 2ds are nonzero 
: 2-N<£<N-2, 2-N<m<N-2, \i - m\ < N - 2. There 
are M = 3N 2 - 9N + 7 such pairs (£,m). Let A dld2d ' A be an M- vector, whose 
components are numbers A d ^ 2ds . Then we have the first key result: A satisfy the 
system of equations (d = d\ + d 2 + d 3 ) 

V 



By moment equations we have created a system of M + d + 1 equations in M 
unknowns. It has rank M and we can obtain unique solution by combination of 
LU decomposition and QR algorithm. The second key result gives us the 2-term 
connection coefficients: 

AA d ld2 = 2 i-d A d ld2 ^ d = dl + d2i A ^ q = a p a q ^ +p 

p 

For nonquadratic case we have analogously additional linear problems for objects 
(10). Also, we use FWT and SSS for computing coefficients of reduced algebraic 
systems. We use for modelling D6,D8,D10 functions and programs RADAU and 
DOPRI for testing. 

As a result we obtained the explicit time solution (7) of our problem. In compar- 
ison with wavelet expansion on the real line which we use now and in calculation 
of Galerkin approximation, Melnikov function approach, etc also we need to use 
periodized wavelet expansion, i.e. wavelet expansion on finite interval. Also in the 
solution of perturbed system we have some problem with variable coefficients. For 
solving last problem we need to consider one more refinement equation for scaling 

JV-1 

function 02 (x): <j>i{x) = J2 a|02(2a; — k) and corresponding wavelet expansion for 

fc=0 

variable coefficients b(t): B° k {b)(f)2{2^x — k), where B 3 k {b) are functionals supported 

k 

in a small neighborhood of 2~^k. 

The solution of the first problem consists in periodizing. In this case we use 

expansion into periodized wavelets defined by (ff^Ax) = 2^ 2 J2 (f>(2^x+2H — k). All 

z 

these modifications lead only to transformations of coefficients of reduced algebraic 
system, but general scheme remains the same. 

II METAPLECTIC WAVELETS 

In this part we continue the application of powerful methods of wavelet analysis 
to polynomial approximations of nonlinear accelerator physics problems. In part 1 
we considered our main example and general approach for constructing wavelet rep- 
resentation for orbital motion in storage rings. But now we need take into account 
the Hamiltonian or symplectic structure related with system (4). Therefore, we 
need to consider generalized wavelets, which allow us to consider the corresponding 
symplectic structures, instead of compactly supported wavelet representation. By 
using the orbit method and constructions from the geometric quantization theory 
we consider the symplectic and Poisson structures associated with Weyl- Heisen- 
berg wavelets by using metaplectic structure and the corresponding polarization. 
In the next part we consider applications to construction of Melnikov functions in 
the theory of homoclinic chaos in perturbed Hamiltonian systems. 



In wavelet analysis the following three concepts are used now: 1). a square inte- 
grable representation U of a group G, 2). coherent states over G, 3). the wavelet 
transform associated to U. 

We have three important particular cases: 
a) the affine (ax + b) group, which yields the usual wavelet analysis 



b). the Weyl-Heisenberg group which leads to the Gabor functions, i.e. coherent 
states associated with windowed Fourier transform. 



In both cases time-frequency plane corresponds to the phase space of group repre- 
sentation. 

c). also, we have the case of bigger group, containing both affine and Weyl-Heisen- 
berg group, which interpolate between affine wavelet analysis and windowed Fourier 
analysis: affine Weyl-Heisenberg group [13]. But usual representation of it is not 
square-integrable and must be modified: restriction of the representation to a suit- 
able quotient space of the group (the associated phase space in that case) restores 

square - integrability: G a wu > homogeneous space. Also, we have more general 

approach which allows to consider wavelets corresponding to more general groups 
and representations [14], [15]. Our goal is applications of these results to problems 
of Hamiltonian dynamics and as consequence we need to take into account symplec- 
tic nature of our dynamical problem. Also, the symplectic and wavelet structures 
must be consistent (this must be resemble the symplectic or Lie-Poisson integrator 
theory). We use the point of view of geometric quantization theory (orbit method) 
instead of harmonic analysis. Because of this we can consider (a) - (c) analogously. 



Let Sp(n) be symplectic group, Mp(n) be its unique two- fold covering - meta- 
plectic group. Let V be a symplectic vector space with symplectic form ( , ), then 
R © V is nilpotent Lie algebra - Heisenberg algebra: 



Sp(V) is a group of automorphisms of Heisenberg algebra. 

Let N be a group with Lie algebra R © V, i.e. Heisenberg group. By Stone- von 
Neumann theorem Heisenberg group has unique irreducible unitary representation 
in which 1 i— > i. This representation is projective: U gi U g , 2 = c(gi,g2) ■ U gi92 , where 
c is a map: Sp(V) x Sp(V) — > S 1 , i.e. c is S^-cocycle. 




P, <P)f\ (x) = exp(i/i(v? -p(x-q))f(x- q) 



Metaplectic Group and Representations. 



[R, V] = 0, [v, w] = (v, w) e R, [V, V] = R. 



But this representation is unitary representation of universal covering, i.e. meta- 
plectic group Mp(V). We give this representation without Stone- von Neumann 
theorem. Consider a new group F — N' tx Mp(V), ix is semidirect product (we 
consider instead of N = R © V the N' = S 1 x V, S 1 = (R/2nZ)). Let V* be dual 
to V, G(V*) be automorphism group of \/*.Then F is subgroup of G(V*), which 
consists of elements, which acts on V* by affine transformations. 
This is the key point! 

Let qi, q n ',Pi, ■■■,p n be symplectic basis in V, a = pdq = J^Pidqi and da be 
symplectic form on V*. Let M be fixed affine polarization, then for a e F the map 
a i— > a gives unitary representation of G: 6 a : H(M) — >■ H(M) 

Explicitly we have for representation of N on H(M): 

(Q q fT(x) = e-^f(x), Q p f{x) = f(x-p) 

The representation of N on H(M) is irreducible. Let A q , A p be infinitesimal opera- 
tors of this representation 

\ = !™ - t [®- tq -I], A p = Jim i[0_ tp - /], 

df 

then A q f(x) = i(qx)f(x), A p f(x) = ^Pj^{x) 

Now we give the representation of infinitesimal basic elements. Lie algebra of the 
group F is the algebra of all (nonhomogeneous) quadratic polynomials of (p,q) 
relatively Poisson bracket (PB). The basis of this algebra consists of elements 
l,9i,-,9n, Pi,-,Pn, QiQj,qiPj, PiPj, i,j = l,...,n, i<j, 

PBis {f,g} = Y / ^^~^^- and {1,0} = /orally 
^ dpj dqi dqi dpi 

{Pi, Qj} = Sij, {piqj, q k } = S ik qj, {piqj,Pk} = SjkPi, {PiPj,Pk} = 0, 
{PiPj, Qk} = kkPj + 5 jk pi, {q^j, q k } = 0, {qiqj,p k } = SikQj ~ SjkQi 
so, we have the representation of basic elements / i— > Af : 1 i— > i, q k i— > ixk, 

This gives the structure of the Poisson manifolds to representation of any (nilpotent) 
algebra or in other words to continuous wavelet transform. 

The Segal-Bargman Representation. 

Let z — l/y/2 ■ (p — iq), z — l/y/2 • (p + iq), p — (p±, ■ ■■,p n ), -F n is the space 
of holomorphic functions of n complex variables with (/, /) < oo, where 



(f,g) = (2n)- n J f(z)g(z)e^ 2 dpdq 



Consider a map U : H — > F n , where H is with real polarization, F n is with complex 
polarization, then we have 

(U*)(a) = J A(a,q)*(q)dq, where A(a, q) = n -^ e -i/2(a 2 W)+V2a q 

i.e. the Bargmann formula produce wavelets. We also have the representation of 
Heisenberg algebra on F n : 

and also : u = dj3 = dp A dq, where (3 = izdz. 



Orbital Theory for Wavelets. 

Let coadjoint action be < g • /, Y >=< f, Ad(gy 1 Y >, where <, > is pairing 
geG, feg*, Y e Q. The orbit is Of = G-f = G/G(f). Also, let A=A(M) be 
algebra of functions, V(M) is A-module of vector fields, A p is A-module of p-forms. 
Vector fields on orbit is 



a(0,X) f (4>) = ±(<i>(exptXf)) 



t=Q 



where G A(0), f G O. Then are homogeneous symplectic manifolds with 
2-form n(a(0, X) f , a{0, Y) f ) =< /, [X, Y] >, and dVt = 0. PB on O have the next 
form {*i,* 2 } = T>i^i)^2 where p is A x (0) -> F(0) with definition fi(p(a),X) 
= i(X)a. Here ^ 2 G and A(0) is Lie algebra with bracket {,}. Now let 

N be a Heisenberg group. Consider adjoint and coadjoint representations in some 
particular case. N = (z,t) G C x R, z = p+iq; compositions in N are (z, t) ■ (V, t') = 
(z + z' ,t + t' + B(z, z')), where B(z,z') =pq' — qp'. Inverse element is (—t,—z). Lie 
algebra n of N is (C, r) G C x R with bracket [((, r), (£', r')] = (0, 5(C, CO)- Centre 
is 5 G n and generated by (0,1); Z is a subgroup exp z. Adjoint representation N on 
n is given by formula Ad(z, t)((, r) = (£, r + £)) Coadjoint: for /en*, # = 
(d ■/)((,() = /(C»t) --Bfo 0/(0,1) then orbits for which /| 5 ^ are plane 
in n* given by equation /(0, 1) = /x . If X = (C, 0), F = (£', 0), X, F G n then 
symplectic structure is 

Q(a(0, X) f , a(G, Y) f ) =< /, [X, F] >= /(0, S(C, C'))/^(C, CO 

Also we have for orbit 0^ = N/Z and is Hamiltonian G-space. 



Kirillov Character Formula or Analogy of Gabor Wavelets. 



Let U denote irreducible unitary representation of N with condition U (0, t) = 
exp(it£) ■ 1, where £ ^ 0,then U is equivalent to representation Te which acts in 
L 2 (R) according to 

Ti(z, t)4>(x) = exp (i£(t + px)) <f>(x — q) 

If instead of N we consider E(2)/R we have S 1 case and we have Gabor functions 
on S 1 . 

Oscillator Group. 

Let O be an oscillator group, i.e. semidirect product of R and Heisenberg group 
N. Let H,P,Q,I be standard basis in Lie algebra o of the group O and H*, P*,Q*, I* 
be dual basis in o*. Let functional f=(a,b,c,d) be al* + bP* + cQ* + dH*. Let us 
consider complex polarizations h = (H, I,P + iQ), h — (I,H,P — iQ) Induced 
from h representation, corresponding to functional f (for a > 0), unitary equivalent 
to the representation 

W(t, n)f(y) = exp(it(h - 1/2)) • U a (n)V(t), 

where V(t) = exp[-it(P 2 + Q 2 )/2a], P = -d/dx, Q = iax, 

and U a (n) is irreducible representation of N, which have the form U a (z) = exp(iaz) 
on the center of N. Here we have: U(n=(x,y,z)) is Schrodinger representation, 
Ut(n) = U(t(n)) is the representation obtained from previous by automorphism 
(time translation) n — > t(n); U t (n) = U(t(n)) is also unitary irreducible repre- 
sentation of N. V(t) = exp(it(P 2 + Q 2 + h — 1/2)) is an operator, which according 
to Stone-von Neumann theorem has the property U t (n) = V {t)U (n)V (t)~ l . 

This is our last private case, but according to our approach we can construct 
by using methods of geometric quantization theory many " symplectic wavelet con- 
structions" with corresponding symplectic or Poisson structure on it. Very useful 
particular spline-wavelet basis with uniform exponential control on stratified and 
nilpotent Lie groups was considered in [15]. 

Ill MELNIKOV FUNCTIONS APPROACH 

In this part we continue the application of the methods of wavelet analysis to 
polynomial approximations of nolinear accelerator physics problems. Now we con- 
sider one problem of nontrivial dynamics related with complicated differential ge- 
ometrical and topological structures of system (4). We give some points of appli- 
cations of wavelet methods from the preceding parts to Melnikov approach in the 
theory of homoclinic chaos in perturbed Hamiltonian systems. 



Routes to Chaos 

Now we give some points of our program of understanding routes to chaos in 
some Hamiltonian systems in the wavelet approach [3]- [11]. All points are: 

1. A model. 

2. A computer zoo. The understanding of the computer zoo. 

3. A naive Melnikov function approach. 

4. A naive wavelet description of (hetero) homoclinic orbits (separatrix) and 
quasiperiodic oscillations. 

5. Symplectic Melnikov function approach. 

6. Splitting of separatrix... — ^stochastic web with magic symmetry, Arnold 
diffusion and all that. 

1. As a model we have two frequencies perturbations of particular case of system 
(4): 

Xi = x 2 x 3 = x 4 , x 5 = 1, x 6 = 1, 

x 2 = — ax i — 6[cos(rx 5 ) + cos(sx 6 )]xi — dx\ — mdx\x\ — px 2 — (p(x 5 ) 
x 4 = ex 3 - f[cos(rx 5 ) + cos(sx 6 )]x 3 - gx\ - kxjx 3 - gx 4 - ip(x 5 ) 

or in Hamiltonian form 

x = J ■ VH(x) + eg(x, 6), 6 = u, (x, 6) e R 4 x T 2 , T 2 = S 1 x S 1 , 

for e = we have: 

x = J- VH(x), & = io (11) 

2. For pictures and details one can see [5], [10]. The key point is the splitting of 
separatrix (homoclinic orbit) and transition to fractal sets on the Poincare sections. 

3. For e = we have homoclinic orbit Xo(t) to the hyperbolic fixed point Xq. For 
e 7^ we have normally hyperbolic invariant torus T e and condition on transver- 
sally intersection of stable and unstable manifolds W s (T e ) and W U (T £ ) in terms of 
Melnikov functions M(O) for Xo(t). 

oo 

M(6) = J VH(x (t)) Ag(x (t),ut + e)dt 

— oo 

This condition has the next form: 



According to the approach of Birkhoff-Smale- Wiggins we determined the region in 
parameter space in which we observe the chaotic behaviour [5], [10]. 

4. If we cannot solve equations (11) explicitly in time, then we use the wavelet 
approach from part 1 for the computations of homoclinic (heteroclinic) loops as 
the wavelet solutions of system (11). For computations of quasiperiodic Melnikov 
functions 

rmT 

M m ^ n (t ) = / DH(x a (t)) A g(x a (t),t + t )dt 
Jo 

we used periodization of wavelet solution from part 1. 

5. We also used symplectic Melnikov function approach 

Mi(z) = lim / {hi,h}w/ tt z)dt 

di(z,e) = h^zl) - hiizl) = eMiiz) + 0(e 2 ) 

where {, } is the Poisson bracket, di(z,e) is the Melnikov distance. So, we need 
symplectic invariant wavelet expressions for Poisson brackets. The computations 
are produced according to part 2. 

6. Some hypothesis about strange symmetry of stochastic web in multi-degree-of 
freedom Hamiltonian systems [11]. 



IV SYMPLECTIC TOPOLOGY AND WAVELETS 

Now we consider another type of wavelet approach which gives us a possibility 
to parametrize Arnold- Weinstein curves or closed loops in Hamiltonian systems by 
generalized refinement equations or Quadratic Mirror Filters equations. 



Wavelet Parametrization in Floer Approach. 

Now we consider the generalization of our wavelet variational approach to the 
symplectic invariant calculation of closed loops in Hamiltonian systems [16]. We 
also have the parametrization of our solution by some reduced algebraical problem 
but in contrast to the general case where the solution is parametrized by construc- 
tion based on scalar refinement equation, in symplectic case we have parametriza- 
tion of the solution by matrix problems - Quadratic Mirror Filters equations [17]. 

The action functional for loops in the phase space is [16] 

F( 7 )= / pdq- f 1 H(t n (t))dt 
Jo 

The critical points of F are those loops 7, which solve the Hamiltonian equations 
associated with the Hamiltonian H and hence are periodic orbits. By the way, all 



critical points of F are the saddle points of infinite Morse index, but surprisingly 
this approach is very effective. This will be demonstrated using several variational 
techniques starting from minimax due to Rabinowitz and ending with Floer ho- 
mology. So, (M, a;) is symplectic manifolds, H : M — * R, H is Hamiltonian, X H is 
unique Hamiltonian vector field defined by 

w(X H (x), v) = -dH(x)(v), v G T X M, x G M, 

where cu is the symplectic structure. A T-periodic solution x(t) of the Hamiltonian 
equations 

x = X H (x) on M 

is a solution, satisfying the boundary conditions x(T) = x(0),T > 0. Let us 
consider the loop space f2 = C°°(S 1 ,R 2n ), where S 1 = R/Zi, of smooth loops in 
R 2n . Let us define a function $ : Q — > R by setting 

rl I rl 

®( x ) = / < ~Ji,x > dt - H(x(t))dt, x G Q 
Jo 2 jo 

The critical points of $ are the periodic solutions of x — Xh{x). Computing the 
derivative at x G Vt in the direction of y G Q, we find 

d f 1 
&(x)(y) = —<S>(x + ey)\ e=0 = J < -Jx - S/H{x),y > dt 

Consequently, <&'(#) (y) = for all y G iff the loop x satisfies the equation 

-Jx(t) - ^H(x(t)) = 0, 

i.e. x(t) is a solution of the Hamiltonian equations, which also satisfies x(0) = x(l), 
i.e. periodic of period 1. Periodic loops may be represented by their Fourier series: 

x (t) = J2 e k2nJt x k , x k G R 2 \ 
kez 

where J is quasicomplex structure. We give relations between quasicomplex struc- 
ture and wavelets in [11]. But now we use the construction [17] for loop parametriza- 
tion. It is based on the theorem about explicit bijection between the Quadratic 
Mirror Filters (QMF) and the whole loop group: LG : S 1 — > G. In particular case 
we have relation between QMF-systems and measurable functions x '■ S 1 — > U(2) 
satisfying 



1 

1 



in the next explicit form 



$i(cj) $i(o; + 7r) 



where 



= xH 



+ 71") 



1 

1 




1 



Also, we have symplectic structure on LG 

1 r 27T 



2, i = 0,l. 



2n Jo 



So, we have the parametrization of periodic orbits (Arnold-Weinstein curves) by 
reduced QMF equations. 

Extended version and related results may be found in [3]- [11]. 
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